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Abstract: 5-aminolaevulinic acid photodynamic therapy (ALA-PDT) is an 
attractive treatment option for nonmelanoma skin tumors, especially for 
multiple lesions and large areas. The efficacy of ALA-PDT is highly 
dependent on the photosensitizer (PS) concentration present in the tumor. 
Thus it is desirable to quantify PS concentration and distribution, preferably 
noninvasively to determine potential outcome. Here we quantified 
protoporphyrin IX (PpIX) distribution induced by topical and intra-tumoral 
(it) administration of the prodrug ALA in basal and squamous cell 
carcinoma murine models by using spatial frequency domain imaging 
(SFDI). The in vivo measurements were validated by analysis of the ex vivo 
extraction of PpIX. The study demonstrates the feasibility of non-invasive 
quantification of PpIX distributions in skin tumors. 
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1. Introduction 

Nonmelanoma skin cancer (basal cell carcinoma (BCC) and squamous cell carcinoma (SCC)) 
is the most frequently diagnosed human cancer. Photodynamic therapy using 5 -aminolevulinic 
acid (ALA-PDT) is an attractive alternative treatment option for nonmelanoma skin cancer in 
situations where surgery is suboptimal, for example due to a large number of tumors or 
tumors located in sensitive areas such as the face. For effective PDT, a sufficient amount of 
photosensitizer (PS) needs to accumulate in tumors. In most clinical studies, PS accumulation 
in each tumor is assumed to be the same for a given administered dose. However, PS 
distribution in tissue can show significant inter and intra-patient heterogeneity [1]. Insufficient 
local PS concentration can lead to suboptimal PDT dose and thus treatment failure. Therefore 
there is considerable motivation for assessing the in vivo PS distribution prior to treatment 
light activation, so that light dose might be adjusted accordingly [2]. 

Fluorescence spectroscopy is frequently used for quantification of the ALA-induced 
protoporphyrin IX (PpIX) content, since the fluorescence contrast is usually higher than the 
absorption contrast in vivo. Spectroscopic measurements, however, do not give information 
about PpIX distribution. Moreover, fluorescence signal is affected by the tissue optical 
properties, and thus is not directly related to PpIX concentration. Ratiometric methods (with 
respect to optical attenuation and/or autofluorescence) partially correct this signal distortion 
[3,4]. Recently spatial frequency domain imaging (SFDI) technique has been proposed for 
quantification of absolute fluorescence concentration by eliminating the variations in 
fluorescence signal due to absorption and scattering at both excitation and emission 
wavelengths [5]. 

In this work, we utilized a custom SFDI system for quantifying the absolute PpIX 
concentration in basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) models, to 
which ALA was applied topically and intra- tumor ally, respectively. PpIX concentration 
mapping clearly showed the improvement in tumor contrast in SCCs, mainly because SCCs 
had high intra-tumor optical absorption compared to surrounding normal tissue. Noninvasive 
imaging results were validated with analysis of the ex vivo extraction of PpIX. Thus we 
conclude that noninvasive SFDI can accurately quantify PpIX concentration distribution 
noninvasively in nonmelanoma skin tumors. 

2. Materials and methods 

2. 1 Animal and tumor models 

All animal experiments followed protocols approved by the IACUC of Roswell Park Cancer 
Institute (RPCI). Transgenic K5-GH mice were acquired from Dr. Andrezj Dlugosz at the 
University of Michigan and bred in our domestic colony. These mice have an activated Sonic 
hedgehog (Shh) signaling pathway and overexpress the Gli2 downstream transcription factor, 
which drives proliferation and leads to the development of spontaneous multiple BCC in the 
skin [6,7]. The spontaneous BCCs are a histologically similar representation of nodular BCCs 
in human skin [6-8]. BCCs can occur with multiple tumors and over large areas, thus topical 
application of ALA (20% ALA in a topical vehicle or cream moisture (Kerastick Dusa, 
Massachusetts)) is the usual clinical option. Two Gli mice from the breeding colony were 
chosen for imaging once they spontaneously developed BCC tumors on their tails. Both mice 
were approximately 6 months old and one mouse had two tumors while the other had one. 
Tumors were selected at the early stage so that tumor sizes were small, ~3 mm in largest 
dimension and reasonably flat. To reduce the curvature effects of mice tails, a piece of clear 
plexiglass was placed over the tail so that the top of the tail and matching Intralipid solution 
surface were kept at the same height. Twenty mg of 5-ALA-hydrochloride (A3785, Sigma) 
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was dissolved in 80 mg of moisturizing cream (Moistrel) immediately before application to 
make a 20% ALA solution. This solution was applied to Gli BCCs 4 h prior to imaging for 
optimal PpIX accumulation and covered with Opsite to keep the cream on the tumor for the 
entire time. On the other hand, squamous cell carcinoma (SCC) can be locally aggressive and 
invasive and may be more suitable for interstitial PDT treatment, with ALA applied 
intravenously (iv) or intratumorally (it). FaDu human epithelial cells derived from a head and 
neck SCC of the hypopharynx were implanted in SCID mice to serve as an SCC model. Two 
SCID mice (age 6 months) were implanted subcutaneously on the lower, right flank with 2.5 x 
10 6 FaDu cells. Flank region was chosen due to its relative flatness as well as to minimize 
breathing artifacts. After 7 days, tumors had grown to ~3 mm diameter in size while 
remaining reasonably flat with negligible thickness. Twenty uL the ALA solution (250 mg/kg 
of ALA, in 0.4 M sodium acetate) were injected directly into the tumor lh prior to imaging. 
Both groups of mice were immobilized during measurements with an inhaled anesthetic 
mixture of isofluorane and oxygen. 

2.2 Custom imaging setup 




Custom lercs tube 
with filters and 
polarizer 



Fig. 1. Diagram and picture of the instrument. The system uses a projector with a built-in red 
LED, 20 nm band-pass filtered at 635 nm light to project images of varying spatial frequencies 
onto tissue simulating phantoms or tissue (Target). The reflected light field is captured by an 
EMCCD camera with adjustable filters to selectively capture emission and excitation light. 
Crossed linear polarizers reduced specular reflection. 

A schematic diagram of the custom spatial frequency domain imaging (SFDI) setup is shown 
in Fig. 1. The setup consisted of a modified projector (Aaxa Tech.) to spatially modulate the 
light into the desired frequency and phase pattern. Light from the projector LEDs was directed 
onto the liquid crystal on silicon (LCOS) chip of the projector where it was spatially 
modulated to produce the desired pattern (Fig. 1). The tube lens of the projector was removed 
and replaced with a 100 mm focal length achromatic lens (#AC254-100-A, Thorlabs) for a 
shorter working distance of 20 cm. The projector contains 3 built-in LEDs (red, green and 
blue). We used two band-pass filters to select the wavelengths of interest from the full 
spectrum of the LEDs. The red LED light of the projector was bandpass-filtered at 630 nm ± 
10 nm for imaging PpIX optical properties at the excitation peak as well as PpIX 
fluorescence. A 660 nm + 20 nm band-pass filter was used to quantify optical properties at the 
emission wavelength of PpIX. The periodic structured illumination pattern with selected 
frequency and phase images was projected onto the target surface. Reflected light was 
collected by the lens and focused to the compact Electron Multiplying CCD (EMCCD) 
camera (Luca, Andor Tech.). Cross-polarizers in front of the projector and camera rejected 
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specularly reflected light, while a 650 nm long-pass fluorescence filter was used for imaging 
PpIX fluorescence and an OD filter for reflectance imaging. 

To quantify absorption and scattering parameters, 6 spatial frequencies ( / ) from 0 to 1.7 

cm" 1 and three phases (0, 27i/3, 4n/3) were used. In reflectance and fluorescence imaging 
modes, EMCCD acquisitions were set at 2s. For optical property quantification, 3 phase and 6 
spatial frequency data acquisition time was ~36s (2s x 6 x 3) for both the excitation and 
emission wavelengths. Fluorescence measurement time was 6s since only the continuous 
wave (CW) frequency (0 cm -1 ) was acquired. The measurement time could be reduced for fast 
acquisition schemes by increasing the LED power (currently it is 100 uW/cm 2 due to the low 
power of the LED inside the projector) and by utilizing the EM gain and removing the OD 



2.3 Model for fluorescence concentration quantification 

PpIX fluorescence concentration was quantified by the Gardner model that corrects raw 
fluorescence signal by compensating for optical absorption (u a ) and scattering (|i s ) loss both 
at excitation and emission wavelengths [5,9]. In this model, the fluorescence correction factor, 
X w {X ex , X em ) is determined by the expression: 



Here X w represents the effective path-length during excitation light penetrating into and 
escape of the emitted fluorescence from the tissue, S = l/[3ju a (ju a +// 5 )] 1/2 is the effective 
penetration depth at the specified wavelength and C l ,C 2 ,C 3 ,k l ,k 2 , k 3 are functions of diffuse 
reflectance R d at excitation and emission wavelengths and are determined empirically 
previously [9]. Bulk optical absorption (ji a ) and scattering (|i s ) parameters were quantified by 
fitting spatial frequency domain reflectance (R d (f)) data with a modified frequency-domain 
diffusion model by using a reference phantom with known optical properties by using custom 
Matlab program [10]. 

2.4 Ex vivo PpIX analysis 

Ex vivo gross quantification of PpIX in tissue samples was performed on solubilized tissue 
extracts [8]. The tumor and neighboring skin pieces were excised and placed in a tube 
containing 0.5 mL of SOLVABLE™ (PerkinElmer). PpIX (Sigma) solution was added to 
control tissue samples solubilized in SOLVABLE™ to obtain a PpIX standard curve. The 
peak fluorescence in the emission spectrum (Ex 407 nm, Em 620-650 nm) was determined as 
a measure of PpIX signal in the solution. Protein content was quantified using the Bio-Rad D c 
protein assay kit and a protein standard curve was created using bovine serum albumin (BSA) 
solutions. PpIX fluorescence signal was normalized to protein levels. 

3. Results and discussion 

3. 1 System calibration 

The system was characterized with multiple phantoms made by titrating the optical properties. 
Similarly, fluorescence phantoms were made by titrating the PpIX concentration. 
Reconstructed raw fluorescence values with respect to concentration were used as a 
"calibration curve" to obtain absolute PpIX concentrations in vivo. 



filter. 




kM ex )/WJ + k 3 (A em )/S(A em ) 
C 2 UJC 3 (A em ) 



(1) 



k 2 {\ x )IS{\ x ) + k 3 {\JI5{\ m ) 
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PplX Concentration (|ig/mL) PplX Concentration (ng/mL) 

Fig. 2. (a) Uncorrected (raw) PplX fluorescence (b) Attenuation corrected PplX fluorescence 
with respect to concentration using all phantoms at different absorption and scattering 
parameters. Inset in b) zoom in to the lowest detected PplX concentration (~4 ng/mL). Error 
bars represent standard deviation of the fluorescence signal due to differences in optical 
properties. 

Figure 2 shows the raw fluorescence intensity with respect to PplX concentrations (0.4, 0.8, 
1.6 ug/mL) at different optical parameter ranges of u a = 0.10 and 0.50 cm" 1 and u s = 10 and 
20 cm" 1 . The raw fluorescence signal showed great variations with respect to optical 
parameters, but attenuation corrected fluorescence showed much less variation with respect to 
PplX concentration. Our phantom experiments showed that the detection sensitivity of the 
system is better than 4 ng/mL by utilization of the electron multiplying gain (set at 50) of the 
EMCCD camera. Higher gain settings and LED power signals can further increase this 
detection sensitivity. 



3.2 In vivo imaging 




Fig. 3. Representative images of a BCC on the tail of a Gli mouse 4h after topical- ALA 
administration, (a) Whitelight structural image showing the tumor area and topical ALA 
application site, (b) Uncorrected PplX fluorescence image showing the tumor and surrounding 
application area, (c) PplX fluorescence concentration indicating higher contrast between the 
tumor and surrounding area compared to the uncorrected image. 



A representative image of a BCC on the tail of a Gli mouse 4 h after topical application of 
ALA on the tumor and surrounding nearby area is shown in Fig. 3. The uncorrected 
fluorescence image (Fig. 3(b)) shows a high PplX signal near the tumor area with some signal 
coming from normal tail tissue (825 + 106 a.u. vs. 465 + 33 a.u. respectively), whereas the 
corrected PplX concentration map shows a more localized signal within the tumor and less 
PplX in the area near tumor (0.44 + 0.07 ug/mL vs. 0.28 ± 0.03 ug/mL respectively). 
Absolute PplX concentrations have a maximum value of -0.5 ug/mL, and an average value of 
-0.35 ug/mL in the region of interest of the tumor. 
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Fig. 4. Representative images of a FaDu tumor lh after it- administration of ALA. (a) 
Whitelight structural image showing the tumor area and ALA injection site, (b) The 
uncorrected fluorescence image does not show localized contrast, (c) PplX fluorescence 
concentration indicating higher contrast between the tumor and surrounding area compared to 
the uncorrected image. 

Figure 4 shows representative images of a FaDu tumor in a SCID mouse lh after it- 
administration of ALA. Figure 4(a) shows the whitelight image to indicate the tumor location. 
Tumors were small (~3 mm diameter), thus the thicknesses were assumed to be negligibly 
small. The uncorrected fluorescence signal (Fig. 4(b)) did not show tumor contrast (1152 + 72 
a.u. vs. 1411 + 86 a.u. in the tumor and normal respectively), but the absolute PplX 
concentration image demonstrated a clear tumor contrast compared to normal surrounding 
tissue (1.10 + 0.12 ug/mL vs. 0.76 + 0.06 ug/mL respectively). As compared to Figs. 3(b) and 
3(c), Figs. 4(b) and 4(c) demonstrate clear improvement in tumor contrast provided by the 
PplX concentration mapping. This is possibly due to the fact that SCCs are well vascularized, 
and have high microvessel density within the tumor compared to BCCs, which have even less 
micro vessel density compared to stroma [11]. Indeed, the contrast in optical properties was 
higher in the FaDu tumor compared to the BCC tumor, as clearly indicated by Table 1. For the 
FaDu tumor shown in Fig. 4, tumor absorption parameter (u a ) was -61% and scattering 
parameter was -12% higher than surrounding normal tissue at the excitation wavelength, and 
-75% and -12% higher at the emission wavelength. This optical contrast was not observed in 
the BCC tumor shown in Fig. 3, where the tumor absorption contrast with respect to skin was 
virtually the same at both the excitation and emission wavelengths, and scattering parameter 
was -8% and -11% lower in the tumor at the excitation and emission wavelengths, 
respectively (Table 1). 

Table 1. Optical Parameters (|ii a , |u s ) Quantified with SFDI at the Emission (-660 nm) and 
Excitation (-630 nm) Wavelengths for Tumor and Normal Sites in SCID and Gli Mice" 









SCID 


Gli 




Wavelength 




FaDu SCC 


normal 


BCC 


normal 


Emission 


Ua 


0.35 ± 0.06 


0.20 + 0.01 


0.21+0.02 


0.22 + 0.02 




Ms' 


9.23 ± 0.99 


8.25 + 0.54 


13.21 + 1.22 


14.62 + 0.70 


Excitation 


Ma 


0.53 ± 0.07 


0.33+0.01 


0.33+0.03 


0.33+0.02 




Ms' 


9.34 + 0.97 


8.32 + 0.50 


13.41+0.63 


14.52 + 0.37 



a Values are given as mean + standard deviation. 



To compare the noninvasive in vivo imaging results with ex vivo PplX analysis, we used a 
hand-drawing tool function (imfreehand, Matlab) to choose not less than 20 x 20 pixel region 
of interest (ROI) for both tumor and normal tissue determined from whitelight and digital 
photographic pictures. Mean and standard deviation values of each ROI were plotted with 
respect to the ex vivo solubilization analysis as shown in Fig. 5. The correlation of ex vivo 
values with respect to uncorrected raw fluorescence signals had a correlation coefficient of 
0.658 (r 2 = 0.658), which improved with the corrected fluorescence signal to r 2 = 0.863. 
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Fig. 5. (a) Uncorrected fluorescence signal vs ex vivo PplX concentration, (b) In vivo PplX 
concentration vs ex vivo PplX concentration. Linearity increased from 0.658 to 0.863. 

Generally, skin layer covers tumor and consideration of skin layer may result in more accurate 
quantification of the tumor tissue. However, in this study, imaged tumors were small and skin 
layer effect is assumed to be minimal compared to human skin. Multi-wavelength approach 
can be helpful to investigate the skin layer, as was demonstrated previously [12]. Since our 
instrument has very few wavelengths, this approach is not practical for accurate 
characterization of skin layers. We are updating our current instrument by adding multi- 
wavelength LED sources and layer effects may be investigated in the future clinical studies. 

4. Conclusion 

We utilized a custom SFDI system to quantify the absolute concentration of PplX in BCC and 
SCC models. PplX concentration quantification showed more contrast improvement for the 
case of SCC model with intra-tumoral administration. Noninvasive imaging results were 
validated by correlation with PplX values determined by the SOLVABLE™ solubilization of 
samples from the same tumors imaged in vivo. Thus we conclude that SFDI can accurately 
quantify PplX concentration distribution noninvasively in nonmelanoma skin tumors. 
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